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Abstract 

We explore the thermodynamic behavior of gases adsorbed within a nanopore. The theoretical 
description employs a simple lattice gas model, with two species of site, expected to describe 
various regimes of adsorption and condensation behavior. The model includes four hypothetical 
phases: a cylindrical shell phase (S*), in which the sites close to the cylindrical wall are occupied, 
an axial phase {A), in which sites along the cylinder's axis are occupied, a full phase {F), in which 
all sites are occupied, and an empty phase {E). We obtain exact results at T = for the phase 
behavior, which is a function of the interactions present in any specific problem. We obtain the 
corresponding results at finite T from mean field theory. Finally, we examine the model's predicted 
phase behavior of some real gases adsorbed in nanopores. 
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I. INTRODUCTION 



A focus of current attention in statistical physics is the behavior of matter in confining 
geometries ^ 0, ^, ||, ||, |^ |, |10|, |lT], |T^ . An extreme version of this problem arises for 
adsorption within nanotubes, a case for which the transverse dimensions may be of order 
molecular sizes. One expects an important parameter in this class of problem to be the ratio 
R* of the diameter of the molecule to that of the tube. When this ratio is of the order one, 
the adsorbate may be well described by a one-dimensional (ID) model. As R* decreases, 
one expects there to arise successively a sequence of onion-like concentric shells of matter; 
the number of possible shells is critically dependent on the value of R*. Accompanying 
the variation in R* is a variation of energy scales, which are the crucial variables in the 
thermodynamics of the system. 

There have been performed many studies of specific geometries and specific adsorbate- 
substrate combinations, as recently reviewed by Gelb et al [|^. However, there have been 
relatively few studies undertaken of the general problem of adsorption in pores in the case 
of variable R*. The present work represents an effort in that direction. Here, we employ 
a highly oversimplified lattice model of adsorption [0, |T5[ designed for cases when one or 
two concentric phases of matter (but no more) may be present. Since the present analysis 
is limited by the assumption of just two distinct species of lattice sites, it describes just the 
R* > 1 regime. Hence, there are assumed to be four possible phases for this geometry: an 
empty phase (E), an axial phase (A), in which atoms are adsorbed only on the cylinder's 
axis, a cylindrical shell phase (S), in which atoms condense close to the cylinder's wall, 
and a full phase (F), in which both axial and shell sites are populated with atoms. These 
are depicted schematically in Fig.l. We assume a model which includes both pore-site 
interactions and nearest neighbor interactions. Since the pore attraction is usually different 
for shell and axial sites, we may think of the axial and shell atoms as two different species 
interacting each other with a common value of the chemical potential fi. The same idea 



was explored in adsorption problems involving two types of binding sites [0. Besides the 
pore attraction, the atoms experience an intra-species interaction (axial-axial or shell-shell) 
and an inter-species interaction (axial-shell). The phase behavior depends on the values of 
these various energies, especially on the attractive or repulsive character of the inter-species 
interaction. 
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Section 2 of this paper presents results at zero temperature (T) for the exact phase 
behavior as a function of the interactions. Section 3 reports a mean field evaluation of the 
phase behavior at finite T. The adsorption behavior given by finite T isotherms is compared 
to the phase diagrams at T = 0. Section 4 describes the relation between the lattice models 
and some examples of possible realistic situations, i.e., gases adsorbed in carbon nanotubes 
of various radii. Ultimately, we would like to relate the systems' properties to energy scales 
present in the real problem. Since these may not be known, it becomes possible in principle 
to deduce these by comparing experimentally observed phase behavior with that predicted 
by the model. In view of the approximations inherent in the lattice model, we believe that 
our results provide a qualitative picture of the expected phase behavior and its evolution 
with the size ratio R* mentioned above. 

II. ZERO TEMPERATURE ANALYSIS 

As a starting point, we consider adsorption in infinite cylindrical pores at T = 0. The 
possible phases are described in the introduction. The cases when the shell-axial interaction 
is attractive or repulsive are investigated separately. We will illustrate in detail our analysis 
for the case of an attractive intcr-species interaction. Initially, to simplify the discussion, 
assume that the analysis can be divided into two alternative approaches. In one, we consider 
the only possible phases to be E, A and F. In the other, we consider just the phases E, S 
and F. We show below that this separation into two distinct treatments encompasses all 
possibilities for the case of attractive interactions between A and S sites. However, in the 
case of a repulsive interaction, this division of the problem into two parts does not work, 
necessitating a somewhat more complicated numerical analysis. 

The phase transition diagrams are constructed on the basis of free energy considerations. 
The shell species is adsorbed on a 2D lattice of sites, with the interaction energy between 
particles at adjacent sites. For simplicity, this 2D lattice on a cylindrical surface is taken as 
a square lattice; hence, the number of nearest shell neighbors of a shell atom Zg is 4. The 
axial species is adsorbed on a ID lattice of sites, of interaction strength and coordination 
number Za = 2. We include also the interaction between axial and shell sites, denoted by 
Sga- Throughout the paper, we express all energies, chemical potentials and temperatures 
in units of e = = , the absolute value of the interatomic interaction. For simplicity. 
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we assume that shell atoms sit on rings whose centers are occupied by axial atoms. The 
number of axial neighbors for a shell site (zsa) is 1 and the number of shell neighbors for an 
axial site (zas) is larger than 1. 

We first determine the equilibrium phase as a function of /x. The axial, shell and full 
grand free energies [Q = F — iJ,N, where F is the Helmholtz free energy) at T = , can be 
written as: 

na^Na{Va-^)-I^N, (l) 
= N,{V, - ^) - llN, (2) 
= NaiVa - |) + N^iVs - |) + N.ZsaSsa - + N^) (3) 

where Na(s) is the number of sites in the axial (shell) phase and Va[s) is the interaction 
potential energy experienced by the axial (shell) site due to the nanotube environment. 
Adsorption in nanopores at T = can occur only if the adsorbate is attracted to the interior 
of the nanopore, i.e. K(s) < 0. We denote the ratio of axial to shell densities (number of 
atoms per pore length) as 7 = Na/Ng. The axial and shell cohesive energies per particle are 
respectively: 

Ea = -(K - Za/2) (4) 

= -{Vs - zj2) (5) 

These energies consist of the pore attraction energy and the nearest neighbor interaction 
(the factor of 1/2 avoids double counting). With this notation, and replacing Zsa by 1, the 
grand free energies can be rewritten as: 

^ = --fEa - (6) 

| = -S.-. (7) 

^ = -7B„- + ,.(7+1) (8) 

One observes that the adsorption behavior (as a function of /i) depends on four parameters: 
7, Esa, Ea and Eg. The T — isotherms are determined by finding the minimum of these Q 
values and comparing the result with the empty lattice result Qe — 0. The axial phase is 
favored relative to the empty phase if fla < 0, i.e. 

/X > -Ea (9) 
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The full phase is lower in grand free energy than the empty phase if < 0, i.e. 



^l> {-Es- Ea-i + esa)/{l + l) 



(10) 



The axial phase is favored relative to the full phase if f^a < implying 



II < -Es + e, 



sa 



(11) 



An analogous argument is true for the shell phase. Jl^ < implies 



At > -E, 



s 



(12) 



Note that < if 



II < -Ea H 



(13) 
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First, we construct two independent phase diagrams with Ea(s) and /i as coordinates, 
corresponding to {E, S, F) and {E, A, F) cases. Then, by inspecting the diagrams, we learn 
how to combine them into a single diagram applicable to both cases at once. We first analyze 
the E, S, F possible phase transitions alone. The n regime of each phase is determined by 
comparison using the equations (10), (12), (13). The transitions between these phases occur 
at values of /i such that the inequalities (10), (12), 13) become equalities. In addition, we 
have to take into consideration that the chemical potential of the pore condensation should 
be smaller than the chemical potential of bulk condensation in the simple cubic lattice Ising 
model, which is — —3. (Of course, transitions can occur within the pore for > //q, but 
one does not ordinarily study them.) Due to this restriction, we can distinguish two cases. 
The first occurs when the S ■<r^ F transition is below saturation {—Ea + egajl < —3). Then, 
all three phases S", F are possible, as shown in Fig.2a. The alternative scenario occurs 
when the S ^ F transition is above saturation {—Ea + esa/7 > ~3). In this case, there are 
only two possible phases, E and S, as shown in Fig.2h. 

The E, A, F phase analysis is very similar to that above for E, S, F. The two cases which 
can be distinguished here are : a) — Eg -\- ega > —3, when all three phases {E, A, F) are 
possible and h) — Es-\- ega < —3 , when there are only two possible phases, E and A. 

So far, the phase transition behavior has been derived from two separate analyses: E, S, F 
and E, A, F. We now show how the parameter values may be assessed in order to establish 
which of the two analyses is appropriate to a given system, i.e. a specified set of parameters. 
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To do so, we need to compare values of and Qg- The difference between the relevant free 
energies satisfies: 

= -Eg + £;„7 - /.(I - 7) (14) 

As can be seen from eqn.9 , the E ■(-^ A transition occurs at /lea — —Ea- In the limit = //go, 
then 

= -{Eg - E^) (15) 

If Ea < Eg, then Qg < Qa at this value of At higher value of fi (> fiea), remains less 
than fla- Hence a transition to the axial phase does not occur for any fi. If, instead. Eg < E^, 
then Q(j < Jl^ and the axial phase is stable at /i — /lea- Is it possible that — Qq changes 
sign for higher /i (corresponding to an A to 5" transition)? This would rec[uire — ^ 
transition value /j, = Has such that 

l^as = — ; (16) 

1-7 

where p — Ea/ Eg < 1. Hence p — 7 > 3(1 — ^)/Ea > 1 — 7. This implies p > 1, which 
violates the assumption Eg < Ea- This rules out such a possibility. 

The same examination can be done at the E S transition line, p^s — ~Eg] we then 
find that for Eg < Ea, the shell phase does not occur. Hence the possibilities are either 
Eg > Ea (never the A phase) or Ea > Eg (never the S phase). This justifies the separate 
analyses used above for the two distinct cases which can arise. 

Because the two cases correspond to different regimes of parameter space. Eg > Ea and 
Eg < Ea, they can be merged in a phase diagram which has as coordinates the interactions 
present in our problem: Ea and Eg. One has only to analyze Fig.2a, b and find the adsorption 
sequences as a function of both interactions when p is increased. Fig.3 exhibits the regimes 
of distinct adsorption sequences. All possible sequences occur except those ruled out by 
the thermodynamic stability condition d/i/dN > 0. The region denoted E corresponds to 
repulsive, or weakly attractive, pore-gas interactions, so that no atoms adsorb inside the 
pore. In the E ^ A region, the shell phase's chemical potential of condensation is greater 
than —3, so the F phase does not occur. Physically, the E ^ A region corresponds to a 
repulsive, or weakly attractive, pore-shell interaction and an attractive pore-axis interaction; 
hence, atoms adsorb only at the axial sites . In the E ^ A ^ F region, the attraction in the 
axial phase is larger than that in the shell phase, so that the axial region is occupied first and 
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then the shell follows at higher /i. Similar reasoning applies to the E ^ S and E ^ S F 
regions. Possibly, the most interesting behavior occurs in the E —>■ F region. In general, as 
seen more clearly at finite T, the axial and shell condensations occur at different chemical 
potentials. However, in the case of an attractive axial-shell interaction, when the shell and 
axial energies per particle are similar, the shell and axial phases become cooperative and 
undergo a common pore filling transition. 

We have examined thus far the case of an attractive axial-shell interaction. In the repul- 
sive case, the inter-species interaction energy (e^a) is positive. Then, we have to take into 
account a new possibility, the transition from axial to shell phase (alone). Physically, this 
means that when the shell atoms are adsorbed, the axial phase, which has a lower density 
than the shell phase, is expelled by the repulsive axial-shell interaction. Therefore, we com- 
pare all the grand free energies Qa, ^s, ^/ with each other and the zero energy of the E phase. 
We present the resulting phase diagrams in {Ea, /i) and {Eg, fi) coordinates in FigAa and 46. 
Both diagrams exhibit all phases and possible transitions E^A^E^S^A^S^A^F 
and S ^ F, but there is no ^ F transition. There are several qualitative differences be- 
tween this case, shown in fig.5, and the attractive interaction case, shown in fig.3. Missing 
in the repulsive case is E ^ F; present in this case are E ^ A ^ S ^ F and E ^ A ^ S 
sequences (absent in the attractive case). The last two are associated with the appearance 
of S, at the expense of A atoms, in order to decrease Q by adding more particles. 



III. FINITE TEMPERATURE ANALYSIS 



In this section, we explore the phase transitions at finite T for a gas within our pore. This 
is a ID system in the thermodynamic limit of divergent length. To study this model, we 
use mean field theory. It is known that ID systems do not exhibit phase transitions at any 
finite T. However, in the present mean field treatment, we obtain a spurious transition. The 



results of an exact calculation of the phase behavior in a square pore were found to be 
qualitatively similar to those of mean field theory, apart from a narrow regime of fi where 
spurious transitions occur in mean field theory; these are replaced by nearly discontinuous 
isotherms in the exact case. We note that gases in some nanoporous media (zeolites or 
nanotube bundles) may represent quasi-lD systems which can go through a genuine phase 
transition when molecules in adjacent pores are coupled. This transition has been studied 



7 



recently in a number of models of gases in pores, by both simulations and exact models ITS 

The occupation probabilities of axial and shell sites are called and n^, respectively. 
We construct the grand free energy of the system and minimize it with respect to Ug and 
Ha- The same procedure was used in Refs.8 and 9 for analyzing layering and wetting phase 
transitions. The energy U of the system is a generalization to finite T of the calculation in 
Section 2. Specifically, the energy is: 

U = NsUsi-^Us + Vs) + NaUai-^Ua + K) + A^s^^s (^^sa'^aCsa) (17) 

and the entropy is written as: 

S = -Ns[ns Inris + (1 - n^) ln(l - n^)] - NaiUalnUa + (1 - ria) ln(l - n^)] (18) 

The minimization of the grand free energy U — TS — fiN with respect to the occupation 
numbers and yields two coupled equations, as found in reference [8]: 

1 



1 + exp(-/5(/i - V; - ZsCsUs - Zsaesaria)) 



= i (19) 

1 + exp(-/3(/i -Va- ZaeaUa - Zasf^saTis)) 

First, we consider the case where the shell- axial interparticle energy e^^ = 0, so that we 
are left with 2 decoupled Ising problems. It is known that a lattice gas can be regarded 
as a lattice of spins, with the conversion s = 2n — 1, J = — e/4 and the magnetic field 
h = ~ V)/2 — ze/A. One can find the chemical potential of condensation from the 
condition for the magnetic transition {h = 0), and the mean field critical temperature Tc 
in the Ising model: jSczJ = 1, where jSc = [ksTc)^^. In the following, we take Boltzmann's 
constant ks = 1. Thus, = ze/A. Therefore, in the decoupled case, the shell and axial 
critical temperatures are Tcs = Zges/A and Tea = Zaea/A, respectively. For simplicity, we 
again use the same axial and shell intra-species interaction eg = €a = e, and scale the 
temperatures with respect to e . Considering a square shell lattice {zg = 4) and a ID axial 
lattice {za = 2), we obtain T^s = 1 and T^a = 0.5. 

Let us consider the effect of turning on the axial-shell interaction. The mean field results 
are shown in fig. 6 for e^a = 1. The chemical potential of condensation is found by a 
Maxwell (equal-area) construction. For a large difference between the energies (per particle) 
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Ea = —{Va — Za/2) and Es = — (K — -2s/2), the shell and axial species behave as in the 
decoupled case; two distinct transitions occur and the transition which occurs first (at lower 
/i) corresponds to a lower free energy. However, in the case of similar energies, the two 
species exhibit a common transition. > 1 in this case because the cooperative system 
behaves like a single species of atoms, with a larger coordination number. 

In order to compare our analysis at finite T with that at T = 0, we keep Vs (or Eg) fixed 
and vary Va (or Ea), so that we move on a line parallel to the Ea axis. In the finite T case, 
we watch the resulting evolution of the axial and shell critical transitions. There arises a 
convenient quantity for characterizing this dependence; this is called 6, defined by: 

S = {Es- esaZsa/2) - {Ea - e,,2;,,/2) (20) 

The evolution of these transitions with 5 is shown in Fig. 7 for three different interaction 
strengths. Consider first the strong attractive case (e^a = !)• For small axial energies per 
particle {6 < —4), the shell condensation occurs at a lower value of fi than that associated 
with full condensation. The axial and shell critical temperatures are the same as in the 
decoupled case (0.5 and 1). This corresponds to the E S ^ F region in figS. When 
6 = —4, the effect of interaction between species becomes significant and the two transitions 
merge. As \6\ approaches 0, the common transition's critical temperature increases to the 
value 1.45 (an increase of 45% ) at 5 = 0. When 6 increases from zero to 4, Tc decreases 
symmetrically with the case 6 < 0. This corresponds to the E ^ F region of {Fig. 3). 
A similar critical temperature dependence on the difference between site binding energies 



was observed in Monte Carlo simulations of benzene condensation in Na-X zeolites HlG 
The difference was in that case, Tc dropped abruptly to zero when 6 exceeded a threshold 
corresponding to a decoupling of the two transitions (since neither species in that case had 
an infinite connected path of its own). When 5 > 4 , the system returns to the case of two 
separate axial and shell transitions. As the axial-shell atractive interaction is reduced, the 
range of 6 values corresponding to cooperative behavior decreases, as shown in Fig.7. Note 
that the maximum value of Tf. for the case e^^ = 0.5 is only 15% greater than that of the 
decoupled shell transition. When ega becomes very small (0.1 in Fig.7), a single transition 
occurs for small \6\, but the transition critical temperature equals that of the shell phase 
alone. 

We have also considered the finite T case of a repulsive interaction, e^^ < 0. Again, we 
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study the behavior with Eg constant and vary Ea, so we move on a hne parallel to the Ea axis 
in Fig.5. The resulting isotherms, corresponding to several different regions in Fig.5, are 
shown in Fig.8. A variety of scenarios can be seen, including those with A either preceding 
or following S. The behavior as a function of Eg is a logical correlate of that shown in Fig.5 
at T = 0. In contrast with the attractive case, there is no E ^ F region, even for similar 
axial and shell energies because the shell atoms, which have a higher density, expel the axial 
atoms. However, there occurs a qualitative similarity of the Tc behavior. At low fi, the axial 
atoms condense first. Then, at higher /i, the shell is occupied while the axis is emptied. 
This transition occurs at the same Tc as the cooperative transition in the attractive case. 
When the external pressure (i.e. /x) is sufficiently high to overcome the axial-shell repulsive 
interaction, a full condensation occurs. These features are expressed in the (T^, /i) diagram 
for various values of e^a and 6 = {Fig. 9). 



IV. REAL GASES IN CARBON NANOTUBES 



We have discussed so far a simple and general theoretical model for adsorption of gases 
in a nanopore. Now, we consider the model's prediction for a specific case - various gases 
adsorbed in C nanotubes. In the spirit of the model, we employ a number of simplifying 
assumptions. The adsorption potential we use is described in ||22| ; it is a sum of Lennard- 
Jones (LJ) two-body interactions between the C atoms (spread into continuous matter) and 
the adsorbate. The energy and distance parameters of this pair potential are obtained from 
semiempirical combining rules involving the LJ parameters of the C atoms (eccO'cc) and 



the adsorbate (€„„, a„„) [El, El, |25 



(^gc = {cTgg + acc)/'2 (21) 
The potential in the nanotube interior at distance r from the axis of the cylinder is PUf : 



V{r,R) = 3n9e,calc[f^{^rM,,{^) - {^fM,{^)] (22) 

where R is the nanotube radius, 9 = 0.32A~^ is the surface density of graphene C atoms 
and 

^"^''^ = /o (l + a;2-2a;cos(0))-/2 ^^^^ 
10 



The adsorption model is simple: the adatoms condense in a close-packed configuration, 
in both the shell and axial phases. We are excluding the case of very large R, which would 
result in the possibility of several concentric shells. As discussed in sections 2 and 3, our 
model has 4 parameters: the shell and axial energies, the ratio of densities (7) and the inter- 
species interaction (e^a). They are not completely independent. One can readily identify 
the axial potential energy as T4 = F(0, R). To find the shell potential, one should examine 
the form of potential. If R is large, V{r, R) has a minimum for a radius Rq larger than 
the hard-core adsorbate radius agg] then it is logical to assume that the gas atoms will be 
adsorbed in the shell phase at this distance {Rg = Rq) and the shell potential is V{Ro, R). 
If the pore radius is small {Rq < Cgg), it is convenient to identify Rg = agg and the shell 
potential Vg = V{agg, R). Geometrical calculations show that this is a good approximation, 
assuming that shell atoms sit near the optimal distance rmin — 2^/^(7gp from axial atoms. 
There is arbitrariness in these assignments, a situation which is inherent in any lattice model. 
Va and Vg lead easily to the axial and shell energies per particle Ea — — (K — Za/2) and 
Es^-{Vs-Zs/2). 

The intra- and inter-species interactions are found using Lennard- Jones parameters for 
the specific gas. The intra-spccics interaction energy is taken as egg and the inter-species 
energy is the adsorbate-adsorbate interaction at r =^ R^ + {Tmin/'^Y- 

esa-^eggii^r-i^r) (24) 

The number of shell atoms contained in a ring of radius Rg is 2'nRs/o-gg and the corre- 
sponding number of axial atoms is 1. Thus, an estimate of the ratio of densities is: 

N a 

7 = ^ = ^ (25) 

Table 1 presents the resulting values of the various parameters for E.2 and Xe inside 
nanotubes of various radii. The sequence of transitions is based on data in Fig.3. We 
note several features of these results. First, the only predicted transition scenarios are 
E^S^F,E^F,E^A and no transition. The E ^ S and E ^ A F sequences 
are not found for H2 or Xe. Physically, E ^ S corresponds to an attractive shell potential 
(negative Vg) but a repulsive axial potential (positive K) ; and E ^ A ^ F corresponds to 
very attractive axial potentials and less attractive shell potentials. These do not occur in 
our model of nanotubes. We do find E ^ A and E ^ S ^ F transitions for a relatively 

11 



large range of nanotube radii. The cooperative behavior E ^ F occurs for a very small 
range of parameters because the gas-gas interaction strength is weak in comparison with the 
nanopore attraction. However, in the case of Xe, which has a much bigger cohesive energy 
{egg = 221K) than H2 {egg = 37K), the mutual transition is more common. The H2 gas 
undergoes the E ^ S ^ F transitions for nanotubes with R > Q A , whereas Xe goes 
through these transitions only for R > 7.7 A. This is due to the difference between these 
molecules' sizes and interaction strengths. For R < 7.3 A, Xe can accommodate only the 
axial phase, whereas the H2 gas would go in the axial phase for R < 5.8 A. For very small 
R (3.5 A for Xe, 3 A for H2), gas does not adsorb at all in nanotubes because the pore gas 
potential becomes repulsive . 

Hartree model calculations and path integral simulations were previously performed for 



adsorption of H2 in C nanotubes of radii 6,7 and 8 A |27|. Our classical results are in 
qualitative agreement with these results. The previous study also found the E ^ S F 
for this range of nanotube radii. However, their quantum calculations allowed them to 
investigate the delocalization of the axial state. For R = 8 A the axial state's probability 
density is no longer confined to the immediate vicinity of the axis, exhibiting a maximum 
near r = 2 A. This is actually not an axial phase, but rather a second shell phase, of small 
radius. In our calculations, the axial phase is confined to the nanotube axis and such a 
second shell phase is not considered. 



V. SUMMARY AND CONCLUSIONS 



We have investigated the adsorption of gases in nanopores, employing a lattice model, 
which we solved exactly at T = and approximately at finite T. Various regimes of transition 
behavior were found, corresponding to a range of interaction strengths. The sequence of 
transitions as a function of fi depends on both the axial-shell interaction energy and the 
difference between the axial and shell energies per particle. When this difference is large, 
the two species condense independently, i.e. the two species are essentially decoupled. When 
this difference is small, the behavior depends on the sign of the axial-shell interaction. For 
Csa > (attractive case), the axial and shell phases undergo a common transition at a 
higher critical temperature. For e^a < 0, an increase of the critical temperature occurs, 
corresponding to an A —> 5* transition. 
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The most important parameter is the radius of the nanotube or, specifically ,the ratio R* 
discussed in the introduction. Even though its value does not appear explicitly in section 2 
and 3, it determines most of the other parameters. This is discussed in section 4, where the 

R dependence of the behavior is explored. Depending on the adsorbate size and interaction 
strength, we find typically that E ^ A occurs for small R*, the E ^ S ^ F occurs for 
large R* and the coupled condensation {E F) occurs for a small range of intermediate 
R*. 

Our approach certainly oversimplifies the real situation in nanopores. First, the lattice 
gas model constrains the atoms to artificial sites that must be identified only by a very 
approximate ansatz, discussed in Section 4. For light gases, such as H2 and He, quan- 
tum effects (such as zero point motion) are very important, yet they are neglected here. 
Nevertheless, we think that our model yields the principal qualitative features of the ad- 
sorption's dependence on the various interactions present in this problem. Thus it should 
help us understand the evolution of adsorption phenomena as a function of adsorbate and 
pore radius. 

This research has been supported in part by grants from the Petroleum Research Fund 
of the American Chemical Society and the Army Research Office. 
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1. Schematic transverse section of a nanotube, showing occupied and unoccupied axial 
and shell sites. 

2. T = phase diagram in the case of an attractive axial-shell interaction, is the 
chemical potential and Eg is defined in Eq.5. Both of these energies are scaled to the intra- 
spccies interaction e. The dashed line is the chemical potential of bulk condensation. We 
distinguish two cases: (a)when the S ■(-^ F transition is present {—Ea + Cso/t > —3) and 
(b)when the S F transition is absent. 

3. T — phase diagram showing the sequence of transitions as a function of shell and 
axial energies, in the case of an attractive axial-shell interaction. Arrows indicate direction 
associated with increasing /i. 

4. T = phase diagram in the case of a repulsive axial-shell interaction: a) as a function 
of Eg and /i, with Ea fixed and b) as a function of Ea and /i, with Eg fixed 

5. T — phase diagram of possible transitions as a function of interactions, in the case 
of a repulsive axial-shell interaction. 

6. Adsorption isotherms in the attractive case: (a) and (c) two transitions Tea = 1, 
Tea = 0.5 (as in the decoupled case) occur for a large difference between axial and shell 
energies; the phase which is first occupied corresponds to a lower energy: axial for (a), shell 
for (c); (b) a cooperative transition at a higher Tc occurs when the axial and shell energies 
are similar. 

7. Dependence of the axial and shell critical temperatures on the difference between 
axial and shell energies. The two transitions for different /i occurring at large 5 merge into 
one common transition when \6\ < Aega- The width of the cooperative behavior regime is 

proportional to e^a- 

8. Isotherms at finite T in the case of repulsive axial-shell interaction: (a) and (c) two 
transitions at the decoupled critical temperatures occur for a large difference between the 
axial and shell energies. The phase which occurs at lower /i corresponds to a lower energy 
per site. In (b) three different transitions occur when the axial and shell energies are similar. 

9. Transition curves in — T plane at 5 = for various values of Cga, in the case of 
a repulsive axial-shell interaction. The axial sites are filled first; then, when the shell gets 
filled, axial atoms are expelled and finally, as increases, the full phase occurs. As in the 
attractive case, the critical temperature is enhanced by the coupling. 
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TABLE I: Possible transitions for different gases and nanotube radii. The Lennard-Jones parame- 
ters are : agg = 3.05 A, egg = 37 K for H2 and agg = 4.1 A, egg = 221 K for Xe. All interaction 
energies {Vs,Va, Eg, Ea, e^n) are expressed in units of the gas hard-core energies Cgg and radii in A. 
The last column shows the sequence of adsorbed phases as /x increases. 
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